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Abstract 

We extend the classical LWR traffic model allowing different maximal 
speeds to different vehicles. Then, we add a uniform bound on the 
traffic speed. The result, presented in this paper, is a new macro- 
scopic model displaying 2 phases, based on a non-smooth 2x2 system 
of conservation laws. This model is compared with other models of 
the same type in the current literature, as well as with a kinetic one. 
Moreover, we establish a rigorous connection between a microscopic 
Follow- The-Leader model based on ordinary differential equations and 
this macroscopic continuum model. 
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1 Introduction 

Several observations of traffic flow result in underlining two different be- 
haviors, sometimes called phases, see pHl [TU [16J, [26] . At low density and 
high speed, the flow appears to be reasonably described by a function of the 
(mean) traffic density. On the contrary, at high density and low speed, flow 
is not a single valued function of the density. This paper presents a model 
providing an explanation to this phenomenon, its two key features being: 

1. At a given density, different drivers may choose different velocities; 

2. There exists a uniform bound on the speed. 

By "bound", here we do not necessarily mean an official speed limit. On 
the contrary, we assume that different drivers may have different speeds at 
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Figure 1: Experimental fundamental diagrams. Left, [281 Figure 1] and, 
right, [261 Figure 1], (see also [23]). 

the same traffic density. Nevertheless, there exists a speed V^nax that no 
driver exceeds. As a result from this postulate, we obtain a fundamental 
diagram very similar to those usually observed, see Figure [Q and Figure [21 
left. Besides, the evolution prescribed by the model so obtained is reasonable 
and coherent with that of other traffic models in the literature. In particular, 
we verify that the minimal requirements stated in [H [13] are satisfied. 

Recall the classical Lighthill-Whitham [32] and Richards [Mj (LWR) 
model 

d tP + d x (pV) = (1.1) 

for the traffic density p. Assume that the speed V is not the same for all 
drivers. More precisely, different drivers differ in their maximal speed w, 
so that V = wijj(p), with w € [w, w],w > 0, being transported along the 
road at the mean traffic speed V. We identify the different behaviors of the 
different drivers by means of their maximal speed, see also [7J [8]. One is 
thus lead to study the equations 

{ d tP + d x (pv)=0 v = w ^ p ), (1.2) 

1 d t w + vd x w = v ; 

Here, the role of the second equation is to let the maximal velocity w be 
propagated with the traffic speed. Indeed, w is a specific feature of every 
single driver, in other words is a Lagrangian marker. Therefore this model 
falls into the class of models introduced in [4J , and later on extended in [3UJ , 
see also [6j formula (1.2)]. 

Introducing a uniform bound Vmax on the speed, we obtain the model 

| dtP + d x (pv) = w . th y = min{y wi){ )] (L3) 

I o t w + v o x w = L J 

We choose to reformulate the above quasilinear system in conservation form, 
similarly to [271 formula (1)], [5J formula (2.2)], [3UJ, formula (1)], see also [36] . 
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as follows: 



| d tP + d ( P v(p, rj)) = w . th = mm f 5 1 

[ d t r) + d x [rj v(p, rj) ) = vr (_ p vr/ J 

(1-4) 

see the Remark l5.3l for further comments on this choice. This model consists 
of a 2 x 2 system of conservation laws with a C 0,1 but not C 1 flow. Note in 
fact that rj/p = w G [w, w]. A 2 x 2 system of conservation laws with a flow 
having a similar C ' 1 regularity is presented in [21] and studied in pp. 

Prom the traffic point of view, we remark that, under mild reasonable 
assumptions on the function ip, the flow in (jl.4[) may vanish if and only 
if p = 0, i.e the road is empty, or p = R, i.e. the road is fully congested. 
It is also worth noting the agreement between experimental fundamental 
diagrams often found in the literature and the one related to (|l,4p . see 
Figure El left. 

From the analytical point of view, we can extend the present treatment to 
the more general case of a maximal speed V max that depends on p, i.e. V max = 
V max (p). However, we prefer to highlight the main features of the model (jl.4p 
in its simplest analytical framework. 

As we already said, the model studied here, inspired from [10J, falls into 
the class of "Aw-Rascle" models. So we could use the approach and the 
theoretical results of [3], which should apply here with minor modifications. 

However, our approach is different: here, contrarily to the above ref- 
erence, we establish directly a connection between the Follow-The-Leader 
model in Section 0] and the macroscopic system (jl.4p , without viewing both 
systems as issued from a same fully discrete system (Godunov scheme) with 
different limits, and without passing in Lagrangian coordinates. For related 
works, including vacuum, see also [2| fTT|. fl8] . 

The present paper is arranged in the following way: in the next section we 
study the Riemann Problem for (jl.4l) and present the qualitative properties 
of this model from the point of view of traffic. In Section [3] we compare 
the present model with others in the current literature and in Section 2] 
we establish the connection with a microscopic Follow-The-Leader model 
based on ordinary differential equations. We also show rigorously that the 
macroscopic model (I1.4D can be viewed as the limit of the microscopic model 
as the number of vehicles increases to infinity. All proofs are gathered in 
the last section. 

2 Notation and Main Results 

We assume throughout the following hypotheses: 
a. R, w, w, Vmax are positive constants, with w < w. 
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b. ip G C 2 ([0,i2]; [0, 1]) is such that 

m = i, vor) 



d 2 

J^(P^(P)) ^ for all p 6 [0, R] . 



C. Zi) > Vrr 



Here, R is the maximal possible density, typically R = 1 if p is normalized 
as in Section HI w, respectively w, is the minimum, respectively maximum, 
of the maximal speeds of each vehicle; V max is the overall uniform upper 
bound on the traffic speed. At [0, the first three assumptions on ^ are 
the classical conditions usually assumed on speed laws, while the fourth one 
is technically necessary in the proof of Theorem 12. 11 The latter condition 
means that all drivers do feel the presence of the speed limit. 
Moreover, we introduce the notation 

F = {(p,w)e[0,R]x[w,w]-.v(p,pw)=V max } (2.1) 
C = {(p,w) E[0,R]x[w,w):v(p,pw) =w^{p)} (2.2) 

to denote the Free and the Congested phases. Note that F and C are closed 
sets and Note also that F is 1-dimensional in the (p, pv) plane 
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Figure 2: The phases F and C in the coordinates, from left to right, (p, pv), 
(p,w) and (p,f]). 

of the fundamental diagram, while it is 2-dimensional in the (p, w) and 
(p, rj) coordinates, see Figure [2j See also Figure [3] to have a vision in three 
dimensions. 

Let be the maximum of the points of maximum of the flow, i.e. = 
max jp € [0, R]: pip(p) = max rg j 0> R] rt/j(r)\. Then, the condition 



wip (p*) > V m 



(2.3) 



is a further reasonable assumption. Indeed, it means that the maximum flow 
is attained in the free phase, coherently with the capacity drop phenomenon, 
see for instance [22]. However, ()2.3p is not necessary in the following results. 

Our next goal is to study the Riemann Problem for (jl.4p . 
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Figure 3: The phases F and C in the coordinates (p,pv,w). Note that F 
is contained in a plane. This figure shows an example of Riemann Problem 
when u l = (p l ,p l v l ,w l ) £ F and u r = (p r , p r v r ,w r ) € C. 



Theorem 2.1 Under the assumptions^, S and\Q, for all states (p\rf), 
{p r ,rf) G F U C, the Riemann problem consisting of with initial data 

ir. \ IP 1 if x < fn , \ rf ifx<0 fn 

p(0, sc) = < % r, _ r/(0, x) = < ' r, _ (2.4) 

rv ; p' «fx > ,v ' y I rf ifx>0 y J 

admits a unique self similar weak solution {p,rj) = (p,rj)(t,x) constructed as 
follows: 

(1) If(p l , V l ),(p r , V r )£F, then 

(2) If (p l ,rf), (p r ,r] r ) G C, then (p,rj) consists of a 1-Lax wave (shock or 
rarefaction) between (p\rf) and {p m ,rf n ), followed by a 2-contact dis- 
continuity between (p m ,rj m ) and (p r ,r] r ). The middle state (p m ,i] m ) is 
in C and is uniquely characterized by the two conditions rf 1 / 1 p m = rf / p l 
and v(p m , rj m ) = v (p r , rj r ). 

(3) If {p l ,rf) G C and (p r ,rf) G F, then the solution {p,rf) consists of a 
rarefaction wave separating (p r ,rf) from a state {p m ,rj m ) and by a linear 
wave separating (p m ,r/ m ) from (p ,t) ). The middle state (p m ,r} m ) is in 
FnC and is uniquely characterized by the two conditions rf 71 / p m = rf / p r 
and v(p m , ry m ) = V. 
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(4) If {p ,rf) € F and (p r \rf) 6 C, then (p,rj) consists of a shock be- 
tween (p l ,rf) and (p m ,ri m ), followed by a contact discontinuity between 
(p m , 7] m ) and (p r , rf). The middle state (p m ,w m ) is in C and is uniquely 
characterized by the two conditions rf 1 / ' p m = rf / p l and v(p m ,rj m ) = 
v{p r ,rf). 

(If (pip(p)} vanishes, then the words "shock" and "rarefaction" above 
have to be understood as "contact discontinuities" ) . 

We now pass from the solution to single Riemann problems to the prop- 
erties of the Riemann Solver, i.e. of the map TZ: (F U C) 2 -» BV(M; C U F) 

such that n((pt,rf), {p r ,r) r )j is the solution to (|1.4j) - (|2.4j) computed at 
time, say, t = 1. 

To this aim, recall the following definition, see [lUj : 

Definition 2.2 A Riemann Solver TZ is consistent if the following two con- 
ditions hold for all {p l ,rf), (p m , rf 1 ), (p r ,r] r ) € F U C , and x € R: 

(CI) IfTz((p l ,r) l ),(p m ,r) m j\ (x) = (p m ,ri m ) and1l({p m ,ri m ),(p r ,T] r ))(x) 
= (p m ,r] m ), then 

(C2) Ifn(ip l , V l ),(p r ,r,n) (x) 

n([ P l ,ri),{p m ^ m )) 
n(( P m , v m ),(p r , v r )) 

Essentially, (CI) states that whenever two solutions to two Riemann prob- 
lems can be placed side by side, then their juxtaposition is again a solution 
to a Riemann problem. Condition (C2) is the vice-versa. 

The next result characterizes the Riemann Solver defined above. 

Proposition 2.3 Let the assumptions^,^ and\Q hold. The Riemann 
Solver TZ defined in Theorem \2.1\ enjoys the following three conditions 

1. R is consistent in the sense of Definition \2.3l 



TZ[{p l ,ri l ),(p m ,ri m U ,if x < x . 
TZ((p m ,rj m ),(p r , V r )) ,ifx>x : 



(p m ,rj m ), then 



K({p l ,n l ),{p r ,rf)) ,ifx<x, 

(p m ,7] m ) , if x > x , 

(p m ,r] m ) , if x < x , 

iz((p l ,rf),(p r ,r, r )) ,ifx>x. 
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Figure 4: The conditions (CI) and (C2). 



2. If (p l ,rf), (p r ,r] r ) £ F, then 1Z y(p l ,rj l ), (p r ,rj r )j is the standard solution 
to the linear system 

\ d t p + 8 X (p Vmax) =0 , , 

\ dtV + d x (r/y max ) = 0, 1 • 1 

3. If (p l ,v l ) € FUC and (p r ,i] r ) G C, then K (V, if), (p r , r) r )) is the 
standard Lax solution to 

d t p + d x (??V(p)) =0 

d tV + d x {^^ P ))=o. (2 - 7) 

Moreover, the conditions ( CI ), 2. and 3. uniquely characterize the Rie- 
mann Solver 1Z. 

The above properties are of use, for instance, in using model (|1.4p on traffic 
networks, according to the techniques described in [15] . 

The next result presents the relevant qualitative properties of the Rie- 
mann Solver defined in Theorem 12.11 from the point of view of traffic. 

Proposition 2.4 Let the assumptions^,^ and\Q hold. Then, the Rie- 
mann Solver TZ enjoys the following properties: 

1. If the initial datum attains values in F, C, or FUC then, respectively, 
the solution attains values in F, C , or F U C . 

2. Traffic density and speed are uniformly bounded. 

3. Traffic speed vanishes if and only if traffic density is maximal. 



4- No wave in the solution to ( f L^[ )- pL^l ) may travel faster than traffic 
speed, i.e. information is carried by vehicles. 
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3 Comparison with Other Macroscopic Models 



This section is devoted to compare the present model (|1 .4[) with a sample 
of models from the literature. In particular, we consider differences in the 
number of free parameters and functions, in the fundamental diagram and 
in the qualitative structures of the solutions. Recall that the evolution 
described by model (jl.4p and the corresponding invariant domain depends on 
the function ip and on the parameters V max _, R, w and w. The fundamental 
diagram of (II. 4p is in Figured! left. 



3.1 The LWR Model 

In the LWR model a suitable speed law has to be selected, analogous 

to the choice of tp in (II. 4p . Besides, in (|1.4jl we also have to set V max , R and 
the two geometric positive parameters w and w. 

The fundamental diagram of (|1.4j) seems to better agree with experi- 
mental data than that of (jl.ip . shown in Figure left. Indeed, compare 




Figure 5: Fundamental diagrams, from left to right, of the (LWR) 
model CEO, of the (AR) model flU]) and of the 2-phase model ([3~2"]) . 

Figure El left with the measurements in Figured! 

As long as the data are in F, the solutions to (jl.4p are essentially the same 
as those of (jl.ip . In the congested phase, the solutions to (|1.4j) obviously 
present a richer structure, for they generically contain 2 waves instead of 1. 
In particular, the (LWR) model (jl.ip may not describe the "homogeneous- 
in-speed" solutions, i.e. a type of synchronized flow, see [261 Section 2.2] 
and [231 [37] , which is described by the 2- waves in (|1.4|) . 

Finally, note that if in (II. 4|) the two geometric parameters w and w coin- 
cide, then we recover the LWR (jl.ip model with V{p) = min{y max , w ip(p)}- 



3.2 The Aw-Rascle Model 

Consider now the Aw-Rascle (AR) model 



dtp + d x 
d t y + d x 



pv(p, y) 
yv{p, y) 



q v(p,y) = ~-p(p) (3.1) 



S 



introduced in |3] and successively refined in several papers, see for instan- 
ce El [HI UHl US ISA [Ml ESI E5] and the references therein. 

Note that w in (jl.4[) plays a role analogous to that of v + p (p) in (13. ip . 

In the (AR) model, i? and the "pressure" function need to be selected, 
similarly to R and ip in (|1.4p . No other parameter appears in (|3.ip . but 
the definition of an invariant domain requires two parameters, with a role 
similar to that of w and w. Indeed, an invariant domain for (13.11) is 



(p, y): p € [0, R] and y € p («_ + , p (u_ + 



see Figure El center, and depends on the speeds t>_ and v+. More recent 
versions of (13. ip contain also a suitable relaxation source term in the right 
hand side of the second equation; in this case one more arbitrary function 
needs to be selected. The original (AR) model does not distinguish between 
a free and a congested phase. However, it was extended to describe two 
different phases in [16]. Further comments on (|3.1|) are found in [29J. 

Concerning the analytical properties of the solutions, the Riemann solver 
for the (AR) model suffers lack of continuous dependence at vacuum, see [8, 
Section 4]. However, existence of solutions attaining also the vacuum state 
was proved in [18], while the 2-phase construction in [16J also displays con- 
tinuous dependence. 

A qualitative difference between the (AR) model and the present one is 
property 3. in Proposition 12.41 Indeed, solutions to (|3.1I) may well have zero 
speed while being at a density strictly lower than the maximal one. 



3.3 The Hyperbolic 2-Phase Model 

Recall the model presented in [TU], with a notation similar to the present 
one: 



Free flow: (p, q) £ F, 
dtp + d x [p- v F (p)] = 0, 
v F (p) = (1 - £) • V 
the phases being defined as 
F = {(p, q) € [0, R] x 



C 



Congested flow: (p, q) £ C, 



I d t p + d x 
| d t q + d x 

vc(p,q) = 



p-vc(p,q)} = 
(q - q*) ■ v c (p,q)] 



■Vf(p) >V f ,q = p- V}, 



(p,q) € [0,R] xR+:v c (p.q) <V C , 



q-q* 



Qi-g* Q2-q* 
R ' R 



(3.2) 



In (|3.2p no function can be selected, on the other hand the evolution depends 
on the parameters V, R and q* while the invariant domains F and C depend 
on Vf, V c , Qi and Q%. A geometric construction of the solutions to (|3.2p in 



the congested phase is in [31 
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The main difference between fundamental diagrams of (13. 2D . see Fig- 
ure right, and that of (jl,4p is that (|3.2p requires the two phases to be 
disconnected: there is a gap between the free and the congested phase. 
This restriction is necessary for the well posedness of the Riemann problem 
for (|3.2p and can be hardly justified on the basis of experimental data. More 
recently, the global well-posedness of the model (|3.2I) was proved in [llj. 

Note that in both models, as well as in that presented in [TB], the free 
phase is one dimensional, while the congested phase is bidimensional. 

The model (|3.2|) allows for the description of wide jams, i.e. of persistent 
waves in the congested phase moving at a speed different from that of traffic. 
Here, as long as ^ (pip(p)) < 0, persistent phenomena can be described 
only through waves of the second family, which move at the mean traffic 
speed. We refer to [29] for further discussions on (|3.2p and comparisons 
with other macroscopic models. 



3.4 A Kinetic Model 

Recall, with a notation adapted to the present case, the kinetic model in- 
troduced in [HJ Section 1]: 



d t r(t, x; w) + d x 



wr(t, x;w)ip\ I r(t,x;w')dw' 



0. (3.3) 



The function ip and the speed w play the same role as here. The unknown 
r = r(t, x; w) is the probability density of vehicles having maximal speed w 
that at time t are at point x. 

In (|3.3p there is one function to be specified, ip, as in (|1.4p . The param- 
eters are R (which is normalized to 1 in [8]), w and w, similarly to (jl.4p . 
Since no limit speed is there defined, no parameter in (|3.3p has the same 
role as here V^ ax . 

Being of a kinetic nature, there is no real equivalent to a fundamental 
diagram for (13. 3p . 

From the analytical point of view, the existence of solutions to (|3.3p 
has not been proved, yet. The main result in [8] only states that (13. 3p can 
be rigorously obtained as the limit of systems of n x n conservation laws 
describing n populations of vehicles, each characterized by their maximal 
speed. 

Let the measure r solve (I3.3P and be such that for suitable functions p 
and w 

r(t,x;-) = p{t,x)S w( t iX ) (3.4) 

where 5 is the usual Dirac measure. Then, formally, (p,w) solves (ll.4p . 
Indeed, for the first equation simply substitute (|3.4|) in (|3.3p and integrate; 
for the second equation substitute (|3.4j) in (|3.3p , multiply by w and integrate 
over \w, w}. 
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Remark that (13. 4p suggests a further interpretation of the quantity p 
in (jl.4p . Indeed, in the present model, at (t, x) vehicles of only one species 
are present, namely those with maximal speed w(t,x). 



4 Connections with a Follow-The-Leader Model 

Within the framework of (jl.3|) , a single driver starting from p at time t = 
follows the particle path p = p(t) that solves the Cauchy problem 

L = v( P (t, Pi t)), w ((t,p(t))) v{ w) = min{ymax) wm (41) 

[p(0)=p 

refer to [12] for the well posedness of the particle path for the LWR model 
(see also [3]). Recall now that w is a specific feature of every single driver, 
i.e. w (t,p(t)) = w(0,p) for all p. On the other hand, from a microscopic 
point of view, if n drivers are distributed along the road, then p is approxi- 
mated by l/{pi + \ — pi), where I is a standard length of a car. 

We fix L > and assume that n drivers are distributed along [—L,L]. 
Then, the natural microscopic counterpart to (II. 4|) is therefore the Follow- 
The-Leader (FTL) model defined by the Cauchy problem 

p* = w (s^> w * = 1 '-"> n 

< p n+ i = V max (4.2) 
Pi (Q) =p i i = 1, . . . ,ra + 1 

where p\ = —L and p n +i = L — I. Proposition 14.11 shows that (|4.2p admits 
a unique global solution defined for every t > and such that pi + \ — pi > I 
for all t > 0. 



Proposition 4.1 Let a., b. and c. hold. Fix L > 0. For any n £ N, 

TO't/i n > 2, choose initial data pf for i = 1, . . . , n satisfying — pf > I. 
Then, the Cauchy problem {l^-lty admits a unique solution pf = pf(t), for 
i = 1, . . . , n + 1, defined for all t > and satisfying p^ +1 (t) — pf(t) > I for 
all t > and for i = 1, . . . , n . 

The proof is postponed to Section 

Our next aim is to rigorously show that in the limit n — > +oo with 
nl = constant > 0, the microscopic model in (|4,2p yields the macroscopic 
one in (|1.4|) . Given the position p % of every single vehicle and its maximal 
speed Wi, for i = 1, . . . , n + 1, the macroscopic variables p, u; are given by 

n I n 
p(x) = > y r r (rc) and = > io? v r r (x) . 
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Note that necessarily pf +1 — pf > I. 

On the contrary, given (p,w) € (L 1 nBV)(M; [0, 1] X [w, w]), with suppp, 
suppu; C [— L, L], we reconstruct a microscopic description defining / = 

p{x) dx^j jn and 

Pn+1 = L — I 

( rPi+i ] 

Pj = max < p € [— L, L]: / = I > for i = 1, . . . , n 

= w(p™+) for i = 1, . . . , n + 1 . 

Note that J" R p(x) dx = nl > 0. Now we are able to rigorously show that, as 
the number of vehicles increases to infinity, the microscopic model in (|4.2p 
yields the macroscopic one in (jl.4p . 

Proposition 4.2 Let a., b. and c. hold. Fix T > 0. Choose (p,w) £ 
(L 1 n BV)(M; [0, 1] x [w,w]) with supp/5, suppu) C [—L,L]. Construct the 

initial data for the microscopic model setting I = (j^ p(x) dx^j jn and 
Pn+1 = L — I 

( rh+i ] 

= max <pG [-L, L]: / p(x)dx = l> fori=l,...,n 

Wi = w(Pi+) for i = 1, . . . , n + 1 . 

Let pf(t), for i = 1, . . . , n, be the corresponding solution to h4-2fy - Define 



n I 

pn{%x) = g^^-p^t) (43) 

n 

w n (t,x) = Yw? Xr „ ,„(*)• (4.4) 

// i/tere exists a pair (p, w) E L°°([0, T]; L 1 (M; [0, 1] x [w, wfj such that 
lim (p n ,w n )(t,x) = (p,w)(t,x) p.a.e. 

n— >+oo 

i/ten, i/te pair (p,pw) is a weak solution to with initial datum (p,pw). 

The proof is postponed to Section \5[ 

5 Technical Details 

We first prove an elementary consequence of our assumption |bl 
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Lemma 5.1 Let ip satisfy \M Then, 
3 p G [0, R[ such that 



if) is constant on [0,p], 
ip is strictly decreasing on [p,R]- 



Proof. Call q{p) = pip{p). If ip is strictly monotone, then p = and the 
proof is completed. Otherwise, assume that ip(pi) = ip(p2) = c for suitable 
pi, P2 G ]0, i?] and p\^ p2- Then, byO, for all p G [pi, P2] we have -0(p) = c 
and g(p) = cp. If ^(0) = c, then the proof is completed. Otherwise, note 
that q'(0) = ip(0) > c contradicts the convexity of q. □ 



Corollary 5.2 Lei ip satisfy^ and\Q Then, 

p < min {p G [0, R\: 3w G [w, w] such that (p, w) G C} . 

The proof is immediate and, hence, omitted. 

In the sequel, for the basic definitions concerning the standard theory of 
conservation laws we refer to [9]. 

Proof of Theorem \2.1\ We consider different cases, depending on the phase 
of the data (|2^|) . 

1. (p l ,r, l ),(p r ,ri r )eF. 

In this case, (|1.4|) reduces to the degenerate linear system (|2.6|) so that 
the problem (ll.4|) ~ (|2.4|) is solved by (j2.5|) . Remark, for later use, that the 
characteristic speed is X F = V m3iX . 

2. (p l , V l ),{p r , v r )eC. 

Now, v(p,r]) = r]ip(p)/p. We show that C is invariant with respect to 
the 2x2 system of conservation laws (|2.7p . To this aim, we compute the 
eigenvalues, right eigenvectors and the Lax curves in C: 



Ai(p, v) 

n(p, v) 

VAi • n 
£i(p\Po,Vo) 



rjtp'(p) +v(p,rf) 



-P 

-v 



-1? [pm] 



Vo- 



Po 



Mp,v) 

r 2 {p,v) 
VA 2 • r 2 

^ 2 (p; p ,Vo) 



v(p,rj) 



tip). 

tip) 





Pv(Po,Vo) 

Hp) 



, po < R. 



When p = R, the 2-Lax curve through. rj ) is the segment p = R, 
77 G [ifou, 

Shock and rarefaction curves of the first characteristic family coincide 
by [5J Lemma 2.1], see also Problem 1, Chapter 5]. The second charac- 
teristic field is linearly degenerate. Hence, (|2.7[) is a Temple system and C is 
invariant, since its boundary consists of Lax curves, see |25l Theorem 3.2]. 
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Thus, the solution to (II, 4D is as described in (2) and attains values in C. 

3. (p l ,r/ 1 ) G C, (p r ,rj r ) G F. 

Let p m satisfy ip{p m ) = V max _p r /rf . Note that such p m exists in ]0, 1[ 
by [b] and [c]. , it is unique by Corollary 15.21 Define rf 71 = (p m /p r )rj r and note 
that (p ,ff), (p m ,i] m ) are connected by a 1-rarefaction wave of (12, 7p having 
maximal speed of propagation Ai (p m ,r/ m ) < Vmax- Hence, a linear wave, 
solution to (12. 6h . can be juxtaposed connecting (p m ,7y m ) to (p l ,r) 1 ) and the 
solution to (|1.4|) is as described at (^J. 

4. (pW) G *\ (pW) G C (see Figureg]). 

Note that system (|2.7h can be considered on the whole of F U C. Also 
this set is invariant for (|2.7p . by |25j Theorem 3.2]. Then, in this case, we 
let (p, ry) be the standard Lax solution to (12. 7p . as described at (4). □ 

Proof of Proposition ^. 31 We consider different cases depending on the phase 
of the data (l2~4"jh 

If (p l ,n l ), (p r ,rj r ) G F, then 7£ (^(p l ,r/ 1 ), (p r ,7? r )^ coincides with the Rie- 
mann solver of a linear system, which satisfies (CI). Condition (C2) is 
immediate since no nontrivial middle state is available. 

Similarly, if (p l , rj l ), (p r , r/ r ) £ C, then 1Z (ip l , rj l ), (p r , rf)^ coincides with 
the standard Riemann solver of a 2 x 2 system, which is consistent. The 
consistency of 1Z then follows by the invariance of C, by 2. in the proof of 
Theorem O 

By the same argument, also the case (p , ry ) G F and (p r ,rf) € C is 
proved. Indeed, in (C2), note that the only possible nontrivial middles 
states are in C. 

Finally, if (p l ,rj l ) G C and (p r ,r] r ) G F, then H ((p l ,<r] 1 ), (p r ,rj r )\ takes 
values in FuC and is the juxtaposition of 2 consistent Riemann problems, 
hence (CI) holds. Concerning (C2), note that the the only possible non- 
trivial middles states are in C, and (C2) follows by the consistency of the 
standard Riemann solver for (|2.7p . 

Thus 1. is proved. Assertions 2. and 3. are immediate consequences of 
the construction of Theorem 12.11 

Assume now that 1Z satisfies 2 and 3. Then all Riemann problems with 
data (pW),(p r ,^) G F, (p l ,rj l ) G F, (p r ,rj r ) G C and (p l , rf), (p r , V r ) G C 
are uniquely solved. The solution to Riemann problems with (p, rf) G C 
and (p r ,rj r ) G F is then uniquely constructed through (CI). □ 

Proof of Proposition \2.4\ Consider the different statements separately. 

1. The invariance of F, C and F U C is shown in the proof of Theo- 
rem [2J3 

2. By the invariance of F U C, it is sufficient to observe that on the 
compact set FUC, the density p, respectively the speed v, is uniformly 
bounded by R, respectively V max . 
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3. It is immediate, see for instance Figure El left. 

4. In phase C we have 



Ai(p,?7) = V^'(p) + v(p,r]) < v(p,rj) and X 2 (p, rj) < v(p, rj). 



In the free phase the wave speed is V max = v(p, rj). The only case left is that 
of a phase boundary with left state in F and right state, say (p r ,r] r ), in C. 
Then, the speed A of the phase boundary clearly satisfies A < \i(p r ,r) r ) < 
v(p r ,rj r ). □ 

Proof of Proposition ^. 1\ Note first that the functions p — * v(p, Wj) in (|4.2|) 
are uniformly bounded and Lipschitz continuous for i = 1, . . . , n. We extend 
them to functions with the same properties and defined on [0, +oo[ setting 



Ui(p) 



^max if 

v(p,Wi) if 
if 



p<0 

p e [o, l] 
P >1. 



(5.1) 



We also note that, for i = 1, . . . , n, the composite applications 5 — ► Ui(l/S), 
can be extended to uniformly bounded and Lipschitz continuous functions 
on [0, +oo[. Now we consider the Cauchy problem 



Pi = 1 

Pn+1 ~ 

pm 



n 



--Pi 



i + l Pi 

max 



1, . . . ,n 



l,...,n + l. 



(5.2) 



Note that p™, for i = 1, . . . , n + 1 are defined in Proposition 14.21 and satisfy 
the condition p^ +1 — p™ > I > 0, for every i = 1, . . . , n. 

By the standard ODE theory, there exists a C 1 solution pf defined as 
long as pf +1 — pf > for all i = 1, . . . ,n. We now prove that in fact 
Pi + i{t) —pf(t) > / for every t > 0. To this aim we assume by contradiction 
that there exist positive i and t*, with i < t*, such that p™ +1 {t) — pf(t) = I 
and Pi +1 (t) - p?(t) < I for every t € )i, t*] . Then, 



Pi(t) = Pi® + />(*) da = P Ut) + T ut 
Jt Jt 



K Pi +1 {s)-Pi( s ) 
This yields a contradiction, since for every t G ] i, £*] 

pr+i(i)-p?(*)>p? + i(*)-p?(*) = ^ 

completing the proof. 



d«=P?(t). 



□ 
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Proof of Proposition ^. 2\ Recall first the definition of weak solution to (11.40 : 
for all cp £ C£°, setting v(p,w) = min{y max , wip(p)}, 




P 
pw 



dtp 



pv(p,w) 
pwv(p,w) 



d x (p \dx dt + 



/ 



P 
pw 



ip(0, x) dx = 



and consider the two components separately. 

Below, 0(1) denotes a constant that uniformly bounds from above the 
modulus of (p and all its derivatives up to the second order. Insert first (|4.3fl 
in the above equality and obtain: 

I n := I I (p n d t ip + p n v{p n ,w n )d x ip) dxdt+ I ptp(0,x)dx 



n 



p?+i(*) 



d t (p + v I 



I 



,Wi \d x (p 



dt 



+ / p n (0,x)(p(0,x)dx + / (p-p n (0,x)) ip(0,x)dx 



n „y 



E 



r7o p?+i(t)-p?(t)j Pn t) 
i 



+E 



Pi+l 



(d t ip{t, x) + pf(t)d x ip(t, x)) dx dt 
p(0,x)dx+ / (p - p n (0,x)) (p(0,x) dx. 



1 Pi+i Pi Jp 

Approximating ip(t,x) with ip(t,pf(t)^ for every x in b?(*)>Pi+i(*)L we 
obtain: 



n r T 



r 



E 



iJo P? + i(t)- P ?(t)J pm dt 



-<p(t,p?(t)) dxdt 



+E 



n „ r p 



I 



P?+l(*) 



i^o P?+i(<)-P?(*)ip?(t) 



O(l) {p? + i(t)-p?(t)) dxdt 



=1 Pi+l Pi J Pi 

n r T 



tp(0,x)dx+ / (p-p n (0,x)) <p(0,x)dx 



(■l j n pT 

ij2 j<p{t, P m)dt+Axj2 o(i){ P ? +1 (t)-put))dxdt 
i=i ^° i=i ^° 



i Pi +1 - p n i jpi 



Pi+i 



ip(0,x)dx+ / (p - p n (0,x)) (p(Q, x) dx 



y 



I 



_ [^(0,^-^(0,^)] 
x Pi+i Pi J Pi 

+0(l)l{pl +1 (T)-pUT))+ [ {p-p n (0,x)) <p{Q,x)dx 
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= 0(1) I (2L + V max T) + [ (p- p n (0,x)) V (0,x)dx 

JR 

and both terms in the latter quantity clearly vanish as n — > +00. 

The computations related to the other component are entirely similar, 
since w is constant along any set of the form 

{(t,x)e[0,T}xR:xe[p?(t), P ? +1 (t)[} 

and the proof is completed. □ 



Remark 5.3 System hi. 2\) is not in conservation form. As far as smooth 
solutions are concerned, it is equivalent to infinitely many 2x2 systems 
of conservation laws. Indeed, introduce a strictly monotone function f E 
C 2 ([w, w); ]0, +00Q . Then, elementary computations show that, as long as 
smooth solutions are concerned, system is equivalent to 

\ dtp + d x (pip{p) g{v/p)) = where V = P /H and 

\ dtV + d x (r]tp(p)g(r]/p)) =0 ' " g(f(w))=w 

Clearly, different choices of f yield different weak solutions to h5.S\) . but they 
are all equivalent when written in terms of p and w. 
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